QCD corrections to double J /if) production in e+e annihilation at = 10.6 GeV 
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Next-to-Leading-Order(NLO) QCD corrections to double J/tp production in e'^e~ annihilation 
at = 10.6 GeV are calculated. We find that they greatly decrease the cross section, with a K 
factor (NLO/LO) ranging from —0.31 ~ 0.25 depending on the renormalization scale. Although 
the renormalization scale dependence indicates a large uncertainty, when combined with the NLO 
QCD corrections to J/tp + rjc production, it can explain why the double J/tp production could not 
be found at B factories while the J/tp + rjc production could, despite the fact that cross section of 
the former is larger than that of the latter at LO by a factor of 1.8. 
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Perturbative quantum chromodynamics calculations 
are essential in the effort to describe large momentum 
transfer processes. To apply it to heavy quarkonium 
physics, the nonrelativistic QCD (NRQCD) factoriza- 
tion approach[l| has been introduced. It allows consis- 
tent theoretical predictions to be made and to be im- 
proved systematically in the QCD coupling constant as 
and the heavy-quark relative velocity v. However, the 
J 1^ polarization measurements at Fermilab Tevatron in 
proton-antiproton collisions [2] and J /ij} production in 
B-factories d, 0: Q have shown that the leading order 
(LO) theoretical predictions in NRQCD could not match 
the experimental results. It is still unclear if all these 
discrepancies can be ameliorated by introducing higher 
order corrections within NRQCD framework. 

One of the most interesting topics in heavy quarko- 
nium physics and NRQCD is the double charmonium 
production in e^e^ annihilation at B factories where 
large discrepancy was encountered. The exclusive pro- 
duction cross section of double charmonium in e^e~ 
J/^ric at ^/s = 10.6 GeV measured by Belle [1, |3] is 
(j[J/^ + B'^'=[> 2] = (25.6 ± 2.8 ± 3.4) fb and by 

BABAR d iscr[J/'(/'+r;c]xB''<=[> 2] = (17.6±2.8;2;i) fb, 
where B^" [> 2] denotes the branching fraction for the rjc 
decaying into at least two charged tracks. Meanwhile, 
the LO NRQCD predictions both in QCD coupling con- 
stant as and the charm-quark relative velocity v, given 
by Braaten and Lee [6], Liu, He and Chao ff], and Hagi- 
wara, Kou and Qiao [S] are only about 2.3 ~ 5.5 fb, which 
is an order of magnitude smaller than the experimental 
results. Ma and Si [§] treated the process by employing 
light-cone method, and similar treatment was performed 
by Bondar and Chernyad [l3|, and Bodwin, Kang and 
Lee In Ref. [13, T(4S') ^ J/V' + 7?c was considered 
by Jia, but this resonant decay contribution to J /ij) + rjc 
cross section turns out to be very small. Such a large 
discrepancy between experimental results and theoretical 
predictions imposes a challenge to the understanding of 
charmonium production based on NRQCD. Many stud- 
ies have been performed to resolve the problem. Braaten 
and Lee Q have shown that the relativistic corrections 
would increase the cross section by a factor of about 2, 



which boosts the cross section to 7.4 fb. The NLO QCD 
corrections to this process have been studied by Zhang, 
Gao and Chao [l^, and also by us in a recent paper [lj|. 
The results show large enhancement to the cross section 
with a K factor (the ratio of the NLO cross section to the 
LO one) of about 2 and the large discrepancy is reduced. 
Moreover, the relativistic corrections have been studied 
by Bodwin, Kang, Kim, Lee and Yu flB'l, and also by He, 
Fan and Chao [l6| , which are also significant. When com- 
bined with the NLO QCD corrections, they may resolve 
the large discrepancy. It seems that the large discrepancy 
between the theoretical result and experimental measure- 
ment is resolved by introducing higher order corrections: 
NLO QCD correction and relativistic correction. 

On the other hand, Bodwin, Lee and Braaten [l7j 
showed that the cross section for the process e"'"e~ 
J /ip + J /"ip may be larger than that for J / ip + rjc hy & 
factor of 1.8, in spite of a suppression factor /aj. that 
is associated with the QED and QCD coupHng constants. 
They suggested that a significant part of the discrepancy 
of J/ip + rjc production may be explained by this pro- 
cess. Hagiwara, Kou and Qiao [1| also calculated and 
discussed this process. And light-cone method is used 
in ref. ^21] by V.V. Braguta. In 2004, a new analysis of 
double charmonium production in e''"e~ annihilation was 
performed by Belle 13] based on a 3 times larger data set 
and no evidence for the process e^e^ J /%p + J was 
found. Both the NLO QCD corrections and relativistic 
corrections to e+e" J /ip + rjc give a large K factor 
of about 2. It is obvious that these two types of correc- 
tions to e+e^ J/ip + J/ip should be studied to explain 
the experimental results. In fact, they have been studied 
by Bodwin, Lee and Braaten for the dominant photon- 
fragmentation contribution diagrams [Toj . The results 
show that the cross section is decreased by K factor of 
0.39 and 0.78 for the NLO QCD and relativistic correc- 
tions respectively. A more reliable estimate, 1.69 ± 0.35 
fb, was given by Bodwin, Lee, Braaten and Yu in ref. [20j . 
In this letter, we present a complete NLO QCD calcula- 
tion to this process and the results show that the cross 
section would be much smaller than the rough estimate 
in Ref. [l^ . Therefore it is easy to understand why there 
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FIG. 1: Feynman diagrams for LO. 



was no evidence for the process e+e^ J/ijj + J/V' 
B-factories. 

At leading order in a^, there are 4 Feynman diagrams. 
Two of them are shown in Fig.[Tl while the other two can 
be obtained by reversing the arrows of the electron lines. 
Momenta for the involved particles are labeled as 

e-(pi) + e+(p2) ^ J/V'(P3) + J/Hp4). (1) 



In the nonrelativistic limit, we can use the NRQCD fac- 
torization formalism to obtain the square of the scatter- 
ing amplitude as shown in Eq. ^ by introducing three 
dimensionless kinematic variables 

. ^ (Pl +P2f £ ^ {Pi ~P3? - ^ (pi -Pif 

where ec and rric are the electric charge and mass of the 
charm quark respectively, and /3 = -^/l — 4/J. -Rs(O) is 
the radial wave function at the origin of J/?/'- The ap- 
proximation Mj/^ = 2toc is taken. After the integration 
of phase space, the total cross section is presented in 
Eq. © . 
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Numerical analysis shows that the leading order total 
cross section is identical to that in ref. ^Ji] by choosing 
the same parameters. 

At next to leading order in as, there are no real emis- 
sion process and we need to calculate only virtual cor- 
rections. Dimensional regularization has been adopted 
for isolating the ultraviolet (UV) and infrared(IR) singu- 
larities. A similar renormalization scheme as in ref. |14l | 
is chosen. The renormalization constants of the charm 
quark mass Zm and field Z2 are defined in the on-mass- 
shell (OS) scheme as 
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where 7^ is Euler's constant, Cp = 4/3 and /i is the 
renormalization scale. The calculation is independent of 
the renormalization scheme of the gluon field and the 
QCD gauge coupling constant since there is no gluon in 
the LO diagrams. 

After having fixed the renormalization scheme and 
omitting diagrams which do not contribute, there remain 
36 NLO diagrams in total, including counter-term dia- 
grams. They are divided into 6 groups as shown in Fig. [21 
UV-divergences only appear in diagrams of group (a) and 



(c), which contain triangle diagrams and corresponding 
counter-term diagrams, and cancel inside both groups. 
Diagrams of all the groups contain IR-divergences. More- 
over, since we take the electron mass nie = in the 
calculation, extra IR-divergences appear in box-diagram 
group (/) , pentagon-diagram group (e) and hexagon di- 
agram group (d) . The combination of soft gluon and zero 
electron mass results in 1/ej^ poles, which however can- 
cel when summing all the diagrams together. Similarly, 
there will be ln(s/mg) divergences in above-mentioned 
groups, and all of them will cancel as expected. In ad- 
dition, diagrams of group (a) and (b) that have a vir- 
tual gluon line connecting the quark pair in J/ip lead to 
Coulomb singularity ~ tt^/i;, which can be isolated by in- 
troducing a small relative velocity v — \pc—Pc\ and taken 
into account by the cc wave function renormalization. 

All the six-point scalar integrals which come from 
hexagon-diagram groups (6) and (d) can be reduced to 
several five-point scalar integrals. Most five-point scalar 
integrals can further be reduced to four-point ones, but 
five of them need to be integrated directly, such as 
Eo[^{pi+P2), ^ip2-pi), ^P3,-5P3,0, 0,0, mc,mc] which 
stems from the reduction of six-point scalar integrals. We 
have developed a complete set of methods to calculate 
tensor and scalar integrals with dimensional regulariza- 
tion, which were realized in our Feynman Diagram Calcu- 
lation package (FDC) [23] . This calculation plays a very 
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important role in the establishment of the methods since 
very difficult scalar integrals are met here. After the 
internal check inside our package on the calculation of 
e+e" J/ipJ/ip was passed, QCD corrections to many 
processes such as J/ip — > ggg, 9gj,jjj are calculated 
and the results are in agreement with previous work in 
Refs. [2^, [l^l- A companion paper about the methods is 
in preparation. All the scalar integrals are calculated an- 
alytically by using FDC. After adding contributions from 
all the diagrams together, the IR-divergent terms cancel 
analytically. The amplitude at NLO can be written as 



(5) 



where M'j^j^q is UV, IR, Coulomb finite and is indepen- 
dent of the renormalization scale /i. Thus our result 
depends on the renormalization scale /i only implicitly 
through as. The expressions for the scalar integrations 
and final results are much more complicated than those 
in Ref. ^4|] for e+e" J/tprjc, since not only the variable 
s but also the variable t appear in the scalar integrations 
of this process. Finally the total cross section at NLO 
can be expressed as: 



dt 
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The analytic expression of K{s,t) is too complicated 
to be presented here. We are content to provide nu- 
merical values for K{s). For rric = 1.5 GeV and s = 
(10.6 GeV)V4m2 = 12.484, K{s) = -11.190 is obtained 
and it implies that the NLO QCD correction deceases the 
cross section. 

In the NLO calculation, we should adopt as in two- 
loop formula with number of active quark flavors Uf = A 
and = 0.338 GeV. The value of the wave function 

MS 

at the origin of J/ip can be extracted from the leptonic 
decay widths: 
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FIG. 2; All Feynman diagrams at NLO in six groups. The 
counter-term diagrams of photon-quark vertex are included 
in (a) and (c) where only the corresponding loop diagrams 
are shown. More diagrams can be obtained by reversing the 
arrows of quark lines and/or interchanging the places of ps 
and p4 and/or interchanging the places of and e~. 
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FIG. 3: Differential cross section as function of \x\ where 
X — cos{9). 9 is the angle between the J/ip and the beam, 
and K = 

a I X I 

of NLO to LO, 



d\x\ 



is the ratio of differential cross section 



is set as 1.5 GeV and fi = is taken. 



By choosing the experimental value given in PDG [25| : 
Fee = 5.55 KeV, together with a = 1/137, Mj/^ = 

2me = 3.0 GeV and as = 0.26, \Rs{0)\^ = 0.944 GeV^ 
is obtained and will be used in the following calcula- 
tion. For other value of rric, it should be multiplied 
by (toc/1.5 GeV)2. We stiU take a = 1/137 and the 
numerical results are showed in Table. [J In Ref. p7j . 
there is an estimate of QCD correction factor K = 
[1 — 8Q!s/(37r)]^ = 0.39 for fragmentation-type diagrams. 
This estimate may be questionable, and the more ap- 
propriate one to the NLO accuracy might instead be 
K = 1 — 4: X 8as/(37r) = 0.15. As a matter of fact, when 



all the diagrams are lumped together, our result confirms 
the latter estimate and disfavors the former one. There- 
fore, together with the contributions from all the other 
diagrams, our results are quite different from their rough 
estimate. A similar plot as in Ref. [13] is shown in 
Fig. [3] to present a more detailed result numerically, and 
it shows that the NLO QCD corrections decrease the dif- 
ferential cross section lesser at the peak around = 1 
and make the peak even sharper. A plot of the cross 
section as function of the center-of-mass (CM) energy is 
presented in Fig. 21 It shows that the cross section at 
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TABLE I: Cross sections with different charm quark mass rric 
and renormahzation scale /x, and -^s = 10.6 GeV. 
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FIG. 4: Cross section as function of the CM energy ^fs with 
nic = 1.5 GeV and ^ — '/s/2. 

NLO becomes smaller as the CM energy goes higher. 
In summary, we have calculated the NLO QCD cor- 



rections to double J/-0 production in e+e" annihilation 
at CM energy 10.6 GeV. Dimensional regularization is 
applied to deal with the UV and IR singularities in the 
calculation, and the Coulomb singularity is isolated by a 
small relative velocity v between the c and c in J /ijj and 
absorbed into the cc bound state wave function. Set- 
ting the electron mass to zero brings extra IR singular- 
ities which cancel in the final results. By taking all the 
diagrams into account, a finite result is obtained. Af- 
ter choosing proper physical parameters, we found that 
the NLO QCD correction K factor ranges from -0.314 to 
0.253 for rric < n < \/s/2. It is strongly dependent on 
the renormahzation scale /i and the dependence indicates 
that the uncertainty is quite large. Therefore it is difficult 
to claim a definitive number for the final result. For the 
default choice of charm quark mass rric = 1.5 GeV and 
renormahzation scale fj, = 2toc, the K factor is 0.077 and 
the NLO cross section of e+e~ J /ip + J/V' is 0.57 fb. 
Meanwhile, the K factor is 1.97 and the NLO cross sec- 
tion is 15.68 fb for e"'"e~ J/ip + rjc given in Ref. p^ . 
Therefore, the cross section of J/ip + rjc production is 
about 25 times larger than that of double J ftp produc- 
tion at NLO. It will make the gap even larger when the 
relativistic corrections are included p^. This large factor, 
albeit with a large uncertainty, can explain why double 
J/f/; production could not be observed while J/ip + rjc 
production could at B-factories. 
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